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Abstract 

An exotic spacecraft propulsion technology is described which ex- 
ploits parity violation in weak interactions. Anisotropic neutrino emis- 
sion from a polarized assembly of weakly interacting particles converts 
rest mass directly to spacecraft impulse. 

1 Introduction 

The photon rocket has long been a familiar example of special relativistic 
kinematics 0]. The speed of light as the ultimate rocket exhaust velocity also 
surfaces in discussions of interstellar travel Q. But few regard the photon 
rocket as even remotely practical. This Note outlines a method of spacecraft 
propulsion which resembles the photon rocket closely, and which appears, in 
principle, capable of reduction to practice. 

Parity nonconservation in the weak interaction can produce thrust by 
anisotropic neutrino emission from a polarized assembly of weakly unstable 
particles or nuclei. A conceptual propulsion system, in which decay or nuclear 
capture of muons from the annihilation of bulk antimatter generates thrust, 
illustrates the principles involved. 

This scheme may be considered a variant on propulsion methods based 
upon reaction from radioactive decay products, which have appeared in the 
literature at intervals since their proposal by Goddard and Tsiolkovsky early 
in this century ||, f§. It is distinct from these in that the decay particles are 
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collisionless and electrically neutral, that the anisotropy of parity-violating 
weak decays resolves the otherwise intractable problem of collimation, and 
that the small cross section for neutrino interactions with other matter per- 
mits thrust generation to take place as a volume, instead of as a surface, 
effect. 

2 Physical Principles 

2.1 Parity nonconservation and thrust generation 

Consider a nucleus of spin J polarized along z which captures a light particle, 
accompanied by the emission of a neutrino with energy Q, momentum q = 
Q/c. The capture imparts to the daughter nucleus momentum q directed 
opposite to the neutrino, along with a small recoil energy. As the only 
preferred axis is the parent nuclear spin J, the rate may be written as a 
Legendre polynomial expansion in cos(6 l ) = q ■ J/qJ: 

J(0)=5>fl(cos(0)) (1) 

Because J is an axial vector, while q is a polar vector, their inner product 
is a psuedoscalar, which is odd under parity. In spherical coordinates, the 
effect of the parity operator is to transform 9 to it — 9. Terms which do not 
conserve parity thus appear as odd harmonics in eqn. [TJ. As a result, the 
average impulse per capture along z 

J 4n dQq z I(9) _ lgai Q 
J 4n dQI{9) 3 c a V c U 

does not vanish. Only a nonvanishing ai coefficient allows neutrinos to trans- 
port a net momentum flux. One concludes that an interaction which con- 
serves parity cannot transport a net flux of momentum from a polarized 
assembly of unstable particles, that a decay whose angular intensity contains 
a linear term in cos(9) will violate parity, and that only such an intensity 
pattern can produce thrust. 

2.2 Microscopic basis of thrust generation 

The weak processes considered here involve muons produced by charged pion 
decay following proton-antiproton annihilation. Muons decay by emission 
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of electron and muon neutrinos in a three-body process. The spectrum of 
free positive muon decay at rest gives < q z >= Q/10 for the average 
momentum carried away by neutrinos, or rj + = 1/10. Negative muons also 
undergo nuclear capture, subsequently emitting neutrons in over 95 % of 
captures. Collimated neutrino emission results from muon capture in which 
either the muon or a nonzero spin target nucleus is polarized. 

The conversion efficiency i]_ for negative muon capture by a polarized 
nucleus is calculated following , incorporating polarization effects after |7j . 
It will be assumed that capture is dominated by L = 1 and L = 2 multipoles, 
that recoil corrections 0(v/c) can be neglected and, initially, that the nucleus 
is 100 % polarized and the muon unpolarized. Nuclear matrix elements for 
individual transitions are computed with a modified closure approximation |3j 
in the pure shell model |J. (The net effect of interference terms is estimated 
to be O (recoil) for medium- mass nuclei.) 

In this manner one may estimate for transitions to a specified final 
angular momentum and parity. Its magnitude typically lies between 0.25 and 
0.3 at the upper end, but can take on values near zero. 

Detailed calculations of muon capture |l(J find the capture strength con- 
centrated in a small number of final spin/parity states. L = 1 transitions 
to states of high spin appear favored, with 65-75 % of the capture going to 
these states. It appears reasonable to assume target nuclei exist with net 
efficiency t]_ in the range 0.15-0.25, within recoil corrections 0(0.04). 

Hyperfine coupling between the target nucleus and an unpolarized IS 
muon will polarize the muon spin slightly at the expense of the target nu- 
cleus spin ||11|| . For high nuclear spins, the polarization loss is small. Muon 
polarization alters the relative importance of transverse and longitudinal mul- 
tipoles |12] so as to increase r]_ somewhat, by perhaps 25- 33 % of its mag- 
nitude for unpolarized muons in favorable cases. 

An estimate for r]_ of ca. 0.2 is taken in the following. 



3 Propulsion System Concept 

The propulsion system operates, in outline, as follows: annihilation yields, 
after free decay of prompt annihilation products, about 3.16 7r ± and 1.85 7r°'s 



Hl3| . It is assumed that neutral annihilation products are discarded so as to 
contribute negligibly to waste heat production. Muons from pion decay are 
decelerated to low energies. Positive muons are stopped and (rapidly) polar- 
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ized, then decay. A polarized nuclear target stops and captures the negative 
muons, resulting in a neutron flux. The neutrons are moderated and captured 
in an isotopic mix of nuclei which regenerates the initial target isotope. Both 
the absorbers which decelerate the muons and the thrust-generating inter- 
action regions are coupled to a high thermal power, low specific mass waste 
heat rejection system. 



3.1 Conversion efficiency 

A simple estimate of pion energetics following annihilation results from the 
pion multiplicities on the assumption that pions share the available mass- 
energy equally. The resulting mean pion kinetic energy is 236 MeV. Pions 
are assumed to decay in flight. Muon neutrinos from charged pion decay 
carry off about 90 MeV on average. It is assumed that charge separation of 
muons may be achieved with gradient drift, and that they are brought to low 
kinetic energy by suitable absorbers. Because of the spread in annihilation 
product kinetic energy, this process requires a degree of momentum sorting. 
The fraction of proton pair mass that appears as muon rest mass is 

(in an obvious notation), divided equally between charge states. Combining 
the conversion efficiency for free positive muon decay with that for negative 
muon capture, adjusted for losses from residual nuclear excitation, leads to 
an overall, ideal efficiency estimate 

77 = ^[r] + + 0.85- 77_] 0.025. (4) 



3.2 Muon interactions 

Stopped positive muons form muonium with unpaired electrons in most sub- 
stances. It is assumed that dynamical polarization by muonium hyperfine 
coupling in a strong magnetic field polarizes muons rapidly compared to their 
lifetime, whereupon their free decay imparts impulse. The thermal energy 
deposited in the ultimate target medium by a stopping muon should not, 
therefore, exceed the hyperfine splitting. 

Negative muons are made to stop in a highly polarized target of a suitable 
nucleus such as 51 V or 55 Mn, whose most probable muon capture daughters 
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are all stable. In most targets negative muons rapidly decelerate to rest, 
whereupon they form muonic atoms. In medium-mass nuclei (Z > 20 or 
so), they overwhelmingly undergo nuclear capture from the muonic IS state 
instead of free decay MM. 



The target may be gaseous ||15| . Negative muons can deposit significant 
heat in the capture target without compromising nuclear polarization, and 
thus only require coarse momentum sorting. In fact, muon capture probably 
must occur in a confined plasma, because hyperfine coupling between K 
and L-shell atomic electrons in the target and the IS" muon-also subject to 
hyperfine coupling with the target nucleus-poses a risk of catastrophic loss 
of target nuclear polarization. A thermal energy of 0.5 keV, or just less than 



5 x 10 6 K, suffices for most medium-mass nuclei [III 



3.3 Neutron issues 

The daughter nucleus deexcites mainly by neutron emission. It is proposed 
to recover the target nucleus for subsequent reuse by operating a catalytic 
neutron capture chain on daughter nuclei. 

Recovery begins with moderation of the neutrons to thermal energies in 
a low-neutron absorption moderator such as graphite or D2 16 0. A jacket 
composed of lighter isotopes of the daughter nucleus captures the neutrons, 
thus running a neutron capture chain terminating in an isobar of the original 
nucleus that quickly decays back to the parent isotope. Chemical separation 
of the parent isotope completes the chain. 

Proton and alpha channels cause some leakage out of the chain. However, 
the nuclei are not lost to the system. Daughter nuclei from (//_, np) or (//_, a) 
reactions lie on the neutron-rich side of stability and will consequently (3- 
decay back to the stability line. An extra neutron capture or two reinserts 
them into the main chain. Neutrons to supply this loss come from diverting 
a fraction of the negative muon flux to capture in heavy nuclides whose 
composition is otherwise immaterial. An estimate based upon the neutron 
multiplicity per muon capture indicates that a 0.01 penalty to rj, of recoil 
order, suffices to redress a 3 % nucleon deficit per capture. 



3.4 Waste Heat Rejection 

For purposes of discussion it will be assumed that waste heat is rejected 
by liquid droplet radiators. Nordley |16| estimates that a specific mass of 
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10~ 5 kg/W t at an outlet temperature of 900 K represents an attainable ex- 
tension of liquid droplet radiator performance. In the present instance, in 
which no penalty accrues to high operating temperatures, it seems permis- 
sible to conjecture that radiators may use as a working fluid a refractory 
metal at an outlet temperature in the low 2000 K's. If evaporation losses are 
low enough, and the nonradiating mass depends weakly upon temperatures, 
then the specific mass of a liquid drop radiator scales as a function of the 
working fluid inlet and outlet temperatures only JTTj . Estimates for a variety 



of candidate working fluids indicate specific masses of order 10 7 kg/W t can 
probably be attained, but extension to greatly lower values appears difficult. 

3.5 Spacecraft dynamics 

Consider a simple rocket propelled by neutrinos. A change in spacecraft rest 
mass dM yields an impulse rjcdM and discarded energy (1 — rj)c 2 dM. Let 
Nt annihilations in a time t impart impulse to the spacecraft at a rate 

p = 2Nm p cr], (5) 

accompanied by local waste heat power 

W = N5W (6) 

within it, where 



SW 



2m p c 2 —n n ± m M c 2 + 7o(m vr ± — m M )c 2 — 7 n 7r om vr oc 2 -|--E • (7) 



The average Lorentz factor 7 is 2.7. E includes ca. 38 MeV from residual 
nuclear excitation, recoil energy from fi + decay, and positron annihilation. 
Waste heat production per annihilation is about 690 MeV. The propulsion 
system mass is 

M = Wa w + m mod + m abs + m p i ant , (8) 

where aw is the specific mass for waste heat rejection, and the masses of neu- 
tron moderator and catalytic chain absorber (governed by neutron scattering 
and absorption lengths, respectively) are, along with those of other propul- 
sion plant components, assumed to scale weakly with power. The high waste 
heat power levels expected in this concept suggest examining the extreme 
in which the propulsion system fixed mass is dominated by the waste heat 
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rejection system. In this limit, the final acceleration is given by the ratio of 
eqns. (W) and (||): 



a = 2 



7]m p c 



erj. 



6Wa w 

With this definition of e the spacecraft expends mass at the constant rate 



(9) 



M _ e/j, f 

Its rapidity p = tanh _1 (t>/c) obeys the traditional rocket equation [1,2], 

p = -7 ] log(l - ^t). 

c 



(10) 



(11) 



Here 1 — efift/c = p(t) is the inverse mass ratio; its final value is M(tf)/Mo = 
Pf. For aw = 10~ 7 kg/W t ,e = 9.1cm/ s 2 . The distance traveled in proper 
time t is 



r^f , ,,dt c 2 
s c i p(fj,(t))— a/i = 



dp 



ep f Ji 



(3{p) d p; 



r/c 2 1 
e p f 



[p f \og(p f ) - n f + 1] , 



(12) 



(13) 



substituting p for v/c = /3, since we are concerned with small rapidities < 0.1. 
In addition to the classical rocket, two other concepts are examined: 



3.5.1 Ram augmentation 

In this variant the spacecraft carries its antiprotonic fuel, but obtains protons 
from the interstellar medium |18 |. Four-momentum balance in the rocket rest 
frame at t = gives 



dM 



d p 



exp(-p) + 2i] - 1' 



(14) 



which integrates to 
\og(p) 



1 - 2r] 



P + log 



exp(— p) + 1 — 2i] 
2rj 



(15) 



a transcendental equation for p(p). The expression for ram deceleration is 
obtained by changing the sign of p in eqn. (|14]). 
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3.5.2 Pellet stream 



The limit of ram augmentation occurs when a stream of fuel pellets is launched 
in advance of the spacecraft with a tailored profile of pellet speeds and launch 
times, so that the spacecraft ingests each fuel pellet at zero velocity in its 
proper frame ||18||. For the pellet stream, 



log(/i) 



1 



2v^ 



rf 



arcsm 



1 + (2r) - 1) cosh(p) 



cosh(p) + (277 - 1) 



(16) 



Modifications to eqn. (||) for ram-augmented systems are second order 
in v/c. For the augmented variants, the distance traveled is computed by 
numerical integration of eqn. (|12"D but with e replaced by e/2 in order to 
maintain consistency with the simple rocket. 



4 Performance Estimates 

Consider a rendezvous mission with a star at a distance of 10 pc, and with a 
duration held fixed at 1000 years. Table 1 displays mission profiles for hybrids 
of pure neutrino rocket, ram- and pellet-augmented rocket, and photon sail 



acceleration at 0.3 gravities p~8[| . 



There is a critical value of e below which it is impossible to travel 10 pc 
in 1000 years. The ram/ram profile requires e c = 23.6. Other profiles in the 
table have smaller critical values. The photon beam/ram deceleration profile 
evidently has the greatest latitude in this regard; for r\ = 0.025, e c = 2.0. If 
the efficiency is as low as rj = 0.015, e c = 3.2. The mass ratio for this last 
example is 2.64. 



5 Discussion 

The propulsion system concept outlined above delivers semirelativistic ra- 
pidity changes for modest mass ratios, but does so at comparatively low 
efficiency. The inefficiency arises from the small charged pion multiplicity 
from annihilation at rest. Purely neutrino- powered examples have been de- 
scribed to illustrate the principle most directly, and for comparison with the 
classical photon rocket. One would doubtless prefer to treat that portion of 
annihilation energy which cannot be converted into collimated muon neu- 
trinos as something other than waste heat; to generate thrust with it! The 
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neutrino rocket may furnish a topping cycle for some other, presumably more 
efficient, method of annihilation-based propulsion. 
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Table 1. Sample Rendezvous Profiles 



Method 


photon/ram 


photon/rocket pellet /ram 


pellet/rocket 


ram/: 


acc. mass ratio 


N/A 


N/A 


2.35 


1.84 


5.13 


acc. t(yr) 


0.11 


0.10 


290 


179 


339 


acc. distance(pc) 


6.1 x 1(T 4 


4.4 x 10~ 4 


1.63 


0.752 


2.10 


max. v/c 


0.035 


0.030 


0.043 


0.030 


0.041 


coast distance(pc) 


9.06 


7.52 


6.91 


5.30 


7.39 


coast t(yr) 


849 


755 


525 


557 


593 


dec. mass ratio 


1.70 


3.28 


1.86 


3.44 


1.82 


dec. t(yr) 


151 


245 


185 


263 


68 


dec. distance(pc) 


0.94 


3.13 


1.46 


3.95 


0.50 


total mass ratio 


1.70 


3.36 


4.37 


6.32 


9.36 



r/o = 0.025. e = 9.0, except for ram/ram (e = 23.6). Photon sail profiles use 
an acceleration of 0.3 gravities. 
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